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ABSTRACT 



If mass and angular momentum are conserved in the cooling flows associated with luminous, 
slowly rotating elliptical galaxies, the flow cools onto extended, massive disks of rotationally 
supported cold gas in the equatorial plane. As the hot interstellar gas approaches the disks, 
its density and thermal X-ray emission increase, resulting in X-ray images that are considerably 
flattened toward the equatorial plane out to an optical effective radius or beyond. Remarkably, the 
flattening of X-ray images due to rotation is very small or absent at the spatial resolution currently 
available to X-ray observations. This is strong evidence that mass and angular momentum are 
not in fact conserved. In particular, if cooling flows are depleted by localized radiative cooling 
at numerous sites distributed throughout the flows, disks of cooled gas do not form and the 
X-ray images appear nearly circular. Even in this case, however, the spatial distribution of the 
cooled gas, and any young stars that may have formed from this gas, would be decidedly flatter 
than the old stellar population; if the young stars are optically luminous, the Balmer lines they 
contribute to the stellar spectra should be more elliptical than the total stellar light. In principle, 
X-ray images of galactic cooling flows can also be circularized by the turbulent diffusion of angular 
momentum away from the axis of rotation. But the effective viscosity of known processes - stellar 
mass loss, supernovae, cooling site evolution, etc. - is insufficient to circularize the X-ray images 
appreciably. Radial gradients in the interstellar iron abundance in elliptical galaxies similar to 
those observed are unaffected by the expected level of interstellar turbulence since these gradients 
are quickly re-established by Type la supernovae. 
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galaxies: interstellar medium - X-rays: galaxies 



galaxies: evolution - galaxies: cooling flows 



1. INTRODUCTION 

Cooling flows in massive elliptical galaxies are 
expected to rotate since most of the interstellar gas 
within an optical effective radius R e is produced 
by mass loss from the slowly rotating stellar sys- 
tems. If the cooling flow proceeds inward conserv- 
ing angular momentum, the interstellar gas ulti- 
mately flows toward a large disk (Rdisk 2i Re) and 
spins up to the local equatorial circular velocity, 
~ 400 km s" 1 (Kley & Mathews 1995; Brighenti 
& Mathews 1996, subsequently referred to as Pa- 



per I) . As interstellar gas approaches the disk, the 
X-ray images are significantly flattened toward the 
equatorial plane on scales ~ R e ~iR e when viewed 
perpendicular to the axis of rotation. For bright 
elliptical galaxies at distance 17 Mpc in the Virgo 
cluster - NGC 4472, NGC 4636, and NGC 4649 - 
R e ~ 7 kpc ~ 1.4'. However, Einstein X-ray im- 
ages show no evidence of rotational flattening ei- 
ther at the resolution of the IPC (~ 55" ~ 4.5 kpc; 
Trinchieri et al. 1986) or at the much higher res- 
olution or HRI (~ 6" ~ 0.49 kpc; Fabbiano, Kim 
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& Trinchieri 1992). In a scries of papers Buotc 
& Canizares have detected large scale X-ray flat- 
tening in several massive elliptical galaxies which 
they attribute to a flattened dark matter poten- 
tial, but the highly flattened inner isophotes antic- 
ipated from momentum-conserving cooling flows 
were not evident (see Buote & Canizares 1998 for a 
review). More recently, Hanlan & Bregman (1999) 
have found little or no rotationally enhanced X- 
ray ellipticity in ROSAT HRI and PSPC images of 
several bright elliptical galaxies. The implications 
of these observations arc not widely appreciated 
among the community of astronomers interested 
in early type galaxies. 

Massive, X-ray luminous elliptical galaxies have 
two principal sources of hot interstellar gas: (1) 
ejection of stellar envelopes from an evolving pop- 
ulation of old stars and (2) inflow of distant cir- 
cumgalactic gas which accumulates in the galactic 
halo either by tidal exchange or secondary infall 
(Brighcnti & Mathews 1999a). If the circumgalac- 
tic gas has a lower specific angular momentum 
than the stars, the X-ray flattening would be less- 
ened as inflowing circumgalactic gas mixes with 
gas lost from the rotating stellar system within 
several R e . But it is unlikely that X-ray disks can 
be eliminated by this means since most of the gas 
in the inner galaxy is produced by stellar mass 
loss and has not flowed in from the galactic halo 
region. Our recent models for the large ellipti- 
cal galaxy NGC 4472 require that at least 60 - 70 
percent of the hot gas in r < R e originates from 
stellar mass loss. This rate of stellar mass ejection 
is expected from normal stellar evolution and is re- 
quired to explain the radial variation of interstellar 
temperature and metallicity typically observed in 
massive elliptical galaxies (Brighenti & Mathews 
1999a). In addition, it is unlikely that the absence 
of observed X-ray flattening can be understood 
simply by the inflow of halo gas into the optical 
centers of elliptical galaxies since it is plausible 
that the extended halo gas has an even larger spe- 
cific angular momentum than the stars. N-body 
simulations of galaxy formation typically produce 
elliptical galaxies and galaxy groups with increas- 
ing specific angular momentum at large radii (e.g. 
Barnes & Efstathiou 1987; Quinn & Zureck 1988). 
But the intrinsic angular momentum in distant 
halo gas is uncertain; gas acquired by tidal dis- 
ruption is expected to have significant rotation, 



but secondary infalling gas that arrives after most 
mergers have occurred may have less net rotation 
than the stellar system. 

However, there is no doubt that the stellar sys- 
tems in most large elliptical galaxies and the gas 
that they expel are rotating. Although the stel- 
lar systems in luminous elliptical galaxies (with 
Lb <L 3 x 10 10 £b,o) are not flattened by rotation, 
their rotation about the minor axis is not small, 
typically 50 - 100 km s _1 (e.g. Binney, Davies, 
& Illingworth 1990), and the gas they expel must 
spin up further as it moves inward. 

These general considerations indicate that the 
net angular momentum of gas ejected from stars 
must be removed to circularize the X-ray images. 
This can be accomplished either by removing gas 
from the flow or by transporting angular momen- 
tum away from the axis of rotation by viscous in- 
teractions. In this paper we explore the impor- 
tance of these two processes and determine their 
relative influence on the X-ray images. The first 
process, localized dropout cooling of interstellar 
gas in regions of low specific entropy, is also re- 
quired to limit the masses of central black holes 
(or dark stellar nuclei) in elliptical galaxies to their 
observed values. Central black holes would be 
about ten times more massive if the cooling flow 
proceeded all the way to the galactic center with 
constant mass flux (Brighcnti & Mathews 1999b). 
The presence of distributed optical emission lines 
from cooling gas also argues for cooling dropout 
(Mathews & Brighcnti 1999a). However, as we 
show below, the strongest argument for the pres- 
ence of localized mass dropout in the interstellar 
gas may be the circularization of X-ray images 
that results. 

We also explore the possibility that interstellar 
viscosity transports angular momentum outward, 
reducing the rotation and X-ray flattening in the 
inner galaxy. An effective viscosity can arise from 
several natural sources: plasma viscosity, interstel- 
lar "turbulent" viscosity driven by stellar mass loss 
or Type la supernovae, the viscosity implicit in re- 
gions of localized interstellar cooling as they sink 
in the galactic potential, and turbulent viscosity 
that may develop from shear instabilities in the 
flow. On several occasions Nulsen and Fabian have 
suggested that the turbulent viscosity is so large 
that the entire interstellar medium becomes ap- 
proximately spherical (Nulsen, Stewart, & Fabian 
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1984; Nulsen & Fabian 1995), but no detailed cal- 
culations were provided to support this conjecture. 

In the following section we describe rotating 
cooling flows in which angular momentum is lost 
by cooling dropout. Then we study the outward 
transport of angular momentum by turbulent vis- 
cosity and the associated diffusion of interstellar 
metallicity. 

2. COOLING DROPOUT IN ROTAT- 
ING COOLING FLOWS 

In the earliest studies of non-rotating cooling 
flows, the interstellar gas was assumed to flow 
completely to the galactic center where, in the 
region of increasing gas density, it cooled by in- 
tense radiative losses. The notion that the hot 
interstellar gas in cooling flows cools not just at 
the centers of elliptical galaxies but throughout 
the galactic volume was first advanced by Fabian, 
Nulsen & Canizares (1982) and Thomas (1986). 
This distributed cooling was introduced to avoid 
the strong central peaking of the X-ray surface 
brightness T, x that occurred in centrally-cooling 
models but not in the galaxies observed. Our re- 
cent studies indicate that the effect of volumetric 
cooling on T, x is more complicated. Although T, x 
is globally reduced if the gas cools before it has 
flowed through the entire potential of the galaxy, 
this reduction is generally rather modest. How- 
ever, in localized regions of intense cooling dropout 
Yi x actually increases because of the enhanced den- 
sity and emission associated with cooling regions; 
if cooling dropout is restricted to a small but finite 
central region, T, x can deviate even further from 
observed profiles. 

Notwithstanding these contrary results, the ev- 
idence for spatially distributed cooling of interstel- 
lar gas is apparent in other observations. The 
dynamically determined masses of central black 
holes (Magorrian ct al. 1998) are typically much 
less than the masses of interstellar gas that have 
cooled over a Hubble time, typically several 10 10 
Mq in massive elliptical galaxies (Brighenti & 
Mathews 1999b). The cooled mass is also far 
larger than the masses of cold HI or H2 gas ob- 
served in luminous elliptical galaxies (e.g., Breg- 
man, Roberts & Giovanelli 1988; Braine, Hcnkcl 
& Wiklind 1997), indirectly indicating that stars 
have formed in the cooled gas. However, if the 



stars produced from cooled gas within R e are op- 
tically dark (e.g., brown dwarfs), the thinness of 
the fundamental plane is difficult to understand 
(Mathews & Brighenti 1999b). This suggests that 
the young, low-mass stars that form are optically 
luminous, which is also consistent with H/3 in- 
dices observed in the stellar spectra in some lumi- 
nous elliptical galaxies (Gonzalez 1993; Faber et 
al. 1995; Trager 1997; Trager et al. 1998; Math- 
ews & Brighenti 1999c). Another indication of dis- 
tributed interstellar cooling is the extended opti- 
cal line emission observed in in the central regions 
of most or all bright elliptical galaxies (e.g., Mac- 
chetto ct al. 1996). We have argued that this gas 
with temperature ~ 10 4 K is cooled interstellar 
gas ionized by stellar UV (Mathews & Brighenti 
1999a). 

In the absence of spatially distributed mass 
dropout or angular momentum redistribution, 
all the cooled gas in rotating cooling flows is 
ultimately deposited in a large disk (Paper I; 
Brighenti & Mathews 1997). Because the gas 
density and radiative emissivity increase as the 
cooling flow approaches the disk, the X-ray images 
in these cooling flow models are highly flattened 
unless the galaxy is viewed along the axis of ro- 
tation. In the following discussion we show that 
spatially distributed mass dropout can sharply re- 
duce, but not eliminate, the rotational flattening 
of X-ray images. As a result, the approximate 
circular shapes of observed X-ray images may be 
the most convincing evidence for the existence of 
cooling dropout. We begin by briefly discussing 
the gravitational potential and stellar dynamics 
in a simple rotating elliptical galaxy; this is fol- 
lowed with a description of a gas-dynamical model 
including the cooling dropout process. 

2.1. Model of a Rotating Galaxy 

For comparison with our earlier models of ro- 
tating elliptical galaxies and for mathematical sim- 
plicity, we adopt here the same E2 model galaxy 
described in detail in Paper I. The ellipsoidal mass 
distributions for stars and dark matter in this 
galaxy are assumed to have King and pseudo- 
isothermal profiles respectively: 

p*(R,z) = p 4l + (R/R„f + (z/z„f}- 3/2 

p h (R, z) = p oh [l + (R/R ch ) 2 + [z/zchf]- 1 - 
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Both distributions are truncated at the same large 
ellipsoidal boundary defined by R t and z t = Rt(l — 
e 2 ) 1 / 2 . The eccentricity of an E2 galaxy (n = 2) 
is e = [1 - (1 - n/10) 2 ] 1 / 2 = 0.6. Scale pa- 
rameters i? c * and R c h for the E2 elliptical can 
be derived from a standard spherical E0 galaxy 
having the same central densities p * and p h by 

(s) 

expanding the characteristic core radii Rc ' for 
the spherical galaxy in the indirection, R c = 
R c s \l — O.ln) -1 / 3 . The physical properties of the 
fiducial non-rotating E0 galaxy are listed in Table 
1. Values of R e , Lb and cr* have been chosen so 
that the E0 galaxy lies on the fundamental plane. 
The total mass of dark matter Mht is nine times 
that of the stellar component, M» t . 

As described in Paper I, we solve the two- 
dimensional Jeans equation to derive the veloc- 
ity dispersion in the meridional plane, a 2 (R, z) = 
er|, = a 2 . Following Satoh (1980), we decompose 
the mean square azimuthal speed into a random 
dispersion cr 2 and a systematic rotation v^ 2 using 
a single parameter k assumed to be constant, 

a\ = v^ 2 -vf = k 2 a 2 + (1 - k 2 )^ 2 . 

When k = 1 the stellar velocity dispersion is 
isotropic and the ellipsoidal distortion is produced 
entirely by rotation; when k = there is no 
net rotation and the ellipsoidal flattening is pro- 
duced entirely by anisotropic velocity dispersion, 
> cr 2 . For a typical slowly rotating giant ellipti- 
cal galaxy in which most (but not all) of the optical 
ellipticity is due to anisotropic orbits, we choose 
k = 0.5 as in the "E2;0.25" model described in Pa- 
per I. From the Jeans equations we can determine 
a self consistent rotation for the stellar system, 
v^ c f > (R, z), which is required for the gas dynamical 
equations. The total galactic potential $(-R, z), 
used in the gas dynamical equation of motion, can 
be derived from the density distributions above 
(see Paper I). The mean stellar temperature that 
appears in the gas dynamical equation for thermal 
energy is given by 

n(R,z) = >^(2a 2 + a 2 ), 

where p m — 0.62 is the mean molecular weight 
and m p is the proton mass. 



2.2. Cooling Dropout 

The physical origin of spatially localized inter- 
stellar cooling and mass dropout is a hypothetical 
ensemble of entropy fluctuations. For example, 
a small region in the cooling flow with low en- 
tropy (in pressure equilibrium with low T, high 
p) will radiate more efficiently and cool before the 
ambient flow reaches the galactic center or cool- 
ing disk. The entropy irregularities have complex 
origins - supernova explosions, stellar mass loss, 
galactic winds and other disturbances in the early 
universe, inhomogeneous magnetic fields, etc. - so 
it is not currently possible to derive them from first 
principles. However, by exploring the influence of 
a variety of assumed mass dropout profiles on the 
X-ray surface brightness and radial mass distri- 
butions in luminous elliptical galaxies, Brighcnti 
& Mathews (1999b) confirmed that the following 
modification of the equation of continuity gives ad- 
equate results: 



The negative term on the right hand side rep- 
resents the disappearance of gas from the flow by 
cooling dropout. This sink term is proportional to 
the local gas density divided by the dropout time 
tdo = 5m p kT/2p m pA, the time for gas to radia- 
tively cool locally at constant pressure. This kind 
of simple mass loss with q = 1, similar to that used 
by Sarazin & Ashe (1989), was shown by Brighenti 
& Mathews (1999b) to be consistent with optical 
and X-ray observations of massive elliptical galax- 
ies. 

The first term on the right side of the equation 
of continuity above describes the source of new 
interstellar gas, which is proportional to the local 
stellar density and to 

a*(i) =a n [t/(t n -U s )]- 1 - 3 , 

the specific mass loss from evolving stars. Here 
a n = 4.7 x 10~ 20 s^ 1 , t n = 13 Gyrs is the cur- 
rent cosmic time and £* s = 1 Gyr is the time that 
the galactic stars are assumed to have formed (see 
Brighcnti & Mathews 1999a for details). We are 
particularly interested in the flattening of X-ray 
images within R e from the centers of large ellipti- 
cals. In this central region most of the interstellar 
gas is produced by stellar mass loss, so we shall 
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generally ignore additional gas flowing in from the 
galactic halo region beyond the optical galaxy. 

Interstellar gas is heated and enriched by Type 
la supernovae which we assume occur at a rate 
SNu(i) = SNu(t„)(i„/i), where SNu(t„) = 0.03 
SNIa per 100 yrs per 10 10 L BQ (Cappellaro 
et al. 1997). The specific rate of mass in- 
jection into the cooling flow by Type la su- 
pernovae is a sn = v sn m sn /M*t where v sn = 
SNu(i)(Z, B /10 10 )/3.15 x 10 9 is the number of su- 
pernovae per second. The rate that gas is supplied 
to the interstellar medium by normal stellar mass 
loss greatly exceeds the mass input from super- 
novae, i.e., a* » a sn . 

2.3. Rotating Cooling Flows with and 
without Dropout 

Figures la and lb illustrate the X-ray surface 
brightness T, X (R, z) in the ROSAT band (0.1 - 2.4 
keV) for the rotating elliptical at time t n = 13 
Gyrs when viewed along the equatorial plane and 
at 60° inclination. As in Paper I, our 2D logarith- 
mic computational grid extends beyond 100 kpc 
in both R and z, but Figure 1 zooms in to show 
only the central ~15 kpc region that is most rele- 
vant to published X-ray images. Figure la shows 
the extremely flattened X-ray appearance of the 
E2 galaxy when no distributed cooling dropout is 
present. This diagram is similar to Figure 12b 
of Paper I with several small differences due to 
the larger value of t n used in that paper and the 
slightly lower Type la supernova rate that we as- 
sume here. When the same rotating cooling flow 
is observed at 60° inclination, the bright X-ray 
emission associated with the outer disk is seen as 
a separate feature at ~ 5 kpc (Figure lb). The 
current ROSAT band luminosity for this galaxy 
within the optical image is L x = 2.0 x 10 40 ergs 
s _1 , about 5 times less than the non- rotating ver- 
sion of the same galaxy (Paper I). However, this 
L x must be regarded as an underestimate because 
our computational proceedure is inaccurate when 
most of the radiative cooling occurs in the central 
zone (Brighenti & Mathews 1999b). 

Figures lc and Id show the same X-ray views 
of the rotating galaxy at t n = 13 Gyrs but now 
with cooling dropout (q = 1) and with all other 
parameters identical to those in Figures la and 
lb. The surface brightness T, X (R, z) in Figures 
lc and Id includes the additional X-ray emission 



from the cooling regions, but the isophotal shapes 
are influenced very little by this contribution. It 
is immediately clear from Figure lc that radiative 
dropout is a major explanation for the nearly cir- 
cular appearance of observed X-ray images. Ev- 
idently, interstellar gas is cooling from the flow 
before it has progressed very far from its stellar 
origin toward the axis of rotation. However, the 
central X-ray bright core of the image is still no- 
ticeably flattened in Figure lc with an X-ray el- 
lipticity that is slightly greater than that of the 
E2 stellar image. (The optical isophotal contours 
would intersect the z-axis at 0.8 of the i?-axis in- 
tersection.) This excess ellipticity of the X-ray 
image is not greatly diminished if the galaxy is 
viewed at an inclination i = 60° (Figure Id). The 
ROSAT luminosity for the q = 1 model E2 galaxy 
at time 13 Gyrs is L x = 8.8 x 10 40 ergs s" 1 . 

When there is no cooling dropout, q = as in 
Figure la, a large and massive disk of cold gas 
forms in the equatorial plane, rotating at the lo- 
cal circular velocity; such disks are described in 
detail in Paper I and in Brighenti & Mathews 
(1997). However, in the calculation with q = 1 
mass dropout by radiative losses as in Figure lc, 
no disk of cold gas forms in the z = plane at least 
beyond ~ 150 pc, the size of our innermost compu- 
tational grid zone. This failure to form a cold disk 
results from the increased radiative dropout rate 
as the density in the cooling flow increases toward 
the equatorial plane, {dp/dt)d oc qp/tdo oc p 2 . 

The azimuthal velocity distributions v ( p(R, z) in 
the interstellar gas shown in Figures 2a and 2b cor- 
respond to the cooling flows imaged in Figures la 
and lc respectively; interstellar rotation is signifi- 
cantly reduced in the presence of cooling dropout. 
Figure 2a shows the azimuthal velocity distribu- 
tion at 13 Gyrs without radiative dropout. As 
described in Paper I, without dropout mass loss 
there is no loss of angular momentum and the 
slow inward velocity of the cooling flow leads to a 
dramatic spinup in angular velocity. In Figure 2a 
the gas is flowing toward a cool disk that extends 
to <~ 5 kpc. As interstellar gas converges toward 
the disk, its azimuthal velocity approaches the 
local circular velocity, ~ 400 km s _1 , as it cools 
and merges with the rotationally supported cooled 
disk. 

With q — 1 (Figure 2b) the rotational veloci- 
ties are reduced by approximately 30 percent in 
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the region shown within 15 kpc of the center. If 
the entire major axis (-5 to +5 kpc) were ob- 
served, the rotational broadening of X-ray emis- 
sion lines would be ~ 400 km s _1 , similar to the 
thermal width at the virial temperature, but reso- 
lution of such lines would be difficult or impossible 
even with the Chandra Observatory. For purposes 
of comparison with future high-resolution X-ray 
observations, we show in Figure 2c the emission- 
weighted line of sight velocity at every position for 
the q = 1 calculation, 



vi os (R,z) = 



J n e n p A R osATV<t> ; n s rfs 
/ n e n p K RO sATds 



where n s is a unit vector along the line of sight 
represented with the s coordinate. The contours 
in Figure 2b and 2c are rather similar due to the 
strong spatial gradient in electron density, but 
vi os (R, z) is about 25 percent smaller than v$(R, z) 
over most of the region plotted. 

Since it is likely that the cooled gas forms into 
optically luminous stars (Mathews & Brighenti 
1999b), it is of interest to compare in Figure 3 
the local density of cooled (dropout) gas pdo(R, z) 
for the q = 1 solution with that of the old stars, 
p*(R,z). Notice that the density distribution of 
cooled gas is considerably flatter than the E2 dis- 
tribution of the old stellar population; the isoden- 
sity ellipticity of the gas just prior to dropout re- 
flects rotational flattening in the galactic poten- 
tial. The orbits of the young stellar population 
are expected to inhabit the same flattened vol- 
ume as the interstellar dropout from which they 
formed, with stellar density contours similar to 
Pdo{R,z). For an optically luminous population 
of young stars formed from cooled gas, we predict 
from Figure 3 that the H/3 equivalent width in the 
stellar spectrum (produced mostly by young stars) 
should exhibit a decidedly greater ellipticity than 
that of the bulk of the stellar light. In addition, 
since the interstellar gas in the cooling flow spins 
up before cooling, we expect that the young stel- 
lar population formed from mass dropout has a 
somewhat greater systemic rotation than that of 
the old stars. 

For a larger dropout coefficient, q = 4, the X- 
ray images are even more circular and the corre- 
sponding azimuthal velocity v${R, z) is lower than 
that of the q = 1 solution. However, when q = 4 
the apparent gas temperature in the cooling flow 



is (unrealistically) lowered by the larger emission 
contribution from cooling sites. We have shown 
that the observed gas temperature profile in NGC 
4472 cannot be fit if q > 4 (Brighenti & Mathews 
1999b). 

3. VISCOUS DIFFUSION OF ANGU- 
LAR MOMENTUM 

We now consider an alternative possible expla- 
nation of the circularity of X-ray images of el- 
liptical galaxies: transport of angular momentum 
away from the axis of galactic rotation by turbu- 
lent diffusion. 

Flow-induced turbulence can occur when the 
ratio of inertial to viscous terms in the equation 
of motion is large, i.e. when the Reynolds number 
1Z = puL/p is much larger than unity. Here p is 
the viscosity and puL represent the characteristic 
density, velocity, and length scale in the cooling 
flow. In an unmagnetized plasma the "molecular" 
viscosity p is proportional to the plasma mean free 
path A oc T 2 /n which can be very large in the hot 
interstellar gas. However, Faraday depolarization 
studies in the interstellar gas in elliptical galaxies 
(Garrington & Conway 1991) indicate microgauss 
magnetic fields for which the proton Larmor ra- 
dius tl is much smaller than A, greatly increasing 
the effective Reynolds number and promoting the 
likelihood of turbulence. To estimate the Reynolds 
number we replace A with the thermal gyroradius 
tl in the expression for the viscosity, p ~ pvth^Li 
where v t h = (Sfc^T/mj,) 1 / 2 is the thermal veloc- 
ity. The resulting Reynolds number for rotating 
cooling flows, 



^pvthrL 



) 



3x10' 



V<j, 



200 km/sy Vkpcy vio 7 K 



( T 



-1/2 



B 
~p~G 



is enormous, so turbulent flow can be expected. 

The Rayleigh criterion establishes a necessary 
condition to avoid shear turbulence in rotating 
flows: for stability to linear perturbations the spe- 
cific angular momentum must increase from the 
rotation axis. By this criterion alone, rotating 
galactic cooling flows appear to be stable. Ini- 
tially, as the gas first enters the flow, its specific 
angular momentum increases outward; this follows 
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from the flat stellar rotation curves typically ob- 
served in slowly rotating massive elliptical galaxies 
(e.g. van der Marel 1991). We have verified that 
the sense of this gradient in the angular momen- 
tum density is preserved in the final cooling flows 
shown in Figure la. 

However, rotating flows that satisfy the Raylcigh 
criterion may still be unstable to non-infinitesimal 
perturbations that result from turbulent mixing. 
A variety of non-linear perturbations are present 
in elliptical galaxies which can generate interstel- 
lar turbulence: transport and dissipation of mass 
ejected from orbiting stars, Type la supernovae, 
the radial sinking of denser cooling sites, buoy- 
ant and shearing magnetic fields, etc. Indeed, 
interstellar mixing (turbulence) due to stellar and 
supernova sources is likely to generate large mag- 
netic fields in the cooling flow gas by the turbulent 
dynamo process. These fields, which can serve to 
redistribute angular momentum, may become dy- 
namically important in the central regions of the 
flow where the magnetic energy density is fur- 
ther increased by the converging flow (Soker & 
Sarazin 1990; Moss & Shukurov 1996; Mathews 
& Brighenti 1997). However, if interstellar turbu- 
lence is invoked to account for the circularization 
of X-ray images, the associated spatial diffusion 
must not suppress abundance gradients observed 
in the hot interstellar gas (Matsushita 1997). 

In the following discussion we estimate the ef- 
fective viscosity corresponding to various sources 
of interstellar mixing or turbulence. Then we de- 
scribe modifications to the gas dynamical equa- 
tions that incorporate a spatially dependent tur- 
bulent viscosity. This is followed by a descrip- 
tion of rotating cooling flows including both mass 
dropout and turbulent viscosity. 

3.1. Sources of Interstellar Viscosity 

The effective viscosity due to stellar mass loss 
and Type la supernovae can be estimated by con- 
structing the viscosity from its simplest kinemat- 
ical representation, p, « piv, where p is the lo- 
cal gas density, I is the mean free path or coher- 
ence length of the disturbed region, and v is the 
characteristic velocity appropriate to each process. 
To estimate the viscosity ^* generated by stellar 
mass loss, we use a length scale w 2a en = 
4.5 x 10 9 p -1 / 3 based on the radius a en of a re- 
cently ejected stellar envelope in pressure equilib- 



rium (Mathews 1990) and the mean turbulent ve- 
locity v*t = 5.7 x 10 5 r^ c 2 cm s" 1 as estimated 
by Mathews and Brighenti (1997). In evaluat- 
ing v*t we use a stellar dispersion cr* = 351 km 
s" 1 (Table 1) and assume that half of the kinetic 
energy in stellar ejecta supports the turbulence. 
With these dimensional values and taking p « 
1.05 x lO" 25 ^- 8 gm cm" 3 we find p,* « 0.06r^ c 4 
gm cm -1 s" 1 for the effective viscosity generated 
by stellar ejecta. Alternatively, the characteris- 
tic length could be based on the total volume 
of interstellar gas disturbed by the motion of an 
ejected stellar envelope. If the stopping distance 
of the stellar envelope after ejection from an orbit- 
ing star is l st « 1.2 x 10 22 rg p 6 c cm (Mathews 1990), 
then t, » £ vol = (7ra 2 „4t) 1/3 « 3.3 x 10 19 r° p 6 c 
cm. In this case the viscosity is somewhat larger, 
M* ~ pv* t tvol w 1.0 r°, p 6 c gm cm" 1 s -1 , reflecting 
the uncertainty in estimating 

For an equally approximate estimate of the 
viscosity corresponding to supernova remnants, 
we use an eddy size i 8n — 3.6 x 10 11 p 1 / 3 cm 
based on the blast wave calculations of Math- 
ews (1990) and a characteristic mixing velocity 
v sn = 2.0 x lO 6 ^- 2 cm s" 1 from Mathews & 
Brighenti (1997), assuming that half of the super- 
nova energy goes into interstellar turbulence. The 
resulting viscosity for Type la supernova explo- 
sions is p sn w 16. 2r^ p 1 c 4 gm cm -1 s _1 . For both 
/U* and p, 8n we ignore a weak time dependence due 
to variable mass loss from evolving stars and the 
time dependent Type la supernova rate respec- 
tively. However, both \x„ and p sn decline with 
galactic radius in a similar fashion due to the de- 
crease in gas and stellar density, so p sn 3> p* is 
expected at all rk pc from these formal estimates. 

An additional viscosity arises from the cooling 
dropout process. The increased density in a cool- 
ing site relative to the surrounding cooling flow 
causes it to sink more rapidly in the galactic po- 
tential than the background flow. The net effect 
of the ensemble of differentially sinking sites is an 
approximately radial exchange of gas in the flow 
with a corresponding viscosity. When a site of 
density p s > p falls at its terminal velocity u t , 
the drag force on the moving site must balance 
the gravitational force, pu 2 irr 2 s ss (4ir /3)r 3 p s <?, 
where g « 5.2 x 10~ 7 r^ pc cm s~ 2 is an ap- 
proximation to the galactic gravitational accel- 
eration. The radius r s of a cooling site can be 
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estimated using the approximate model for an 
isobaric, steady state cooling site described by 
Mathews & Brighenti (1999a). Within the cool- 
ing site structure, the temperature varies nearly 
linearly with radius, T = \KK.(T)Y' z r where 
KK = (8Tr/5){p m /k) 3 P 2 m p fA(T)/m s = 5.2 x 
10 _38 r fep 3 c 6 . Here p m = 0.62 is the molecular 
weight, f(p m ) = 0.58, and P«1.9x 10- 10 r^ c 8 
dyne is the approximate interstellar pressure vari- 
ation. The flow into each site is m s = 10 -6 M Q 
yr _1 since line emission from ~ 10 6 sites is re- 
quired to match the total H/3 luminosity from lu- 
minous elliptical galaxies (Mathews & Brighenti 
1999a) and approximately 1 M Q of interstellar 
gas cools each year. Taking T « at the site 
center, the radius where T = 10 6 K, at which 
the site density is Wp, is r s w 2.7 x 10 18 r^ p 2 c 
cm and the corresponding terminal velocity is 
ut = [Agr s {p s / P )/i\ 1 / 2 = 43 r kpc 01 km s" 1 . The 
effective viscosity is then p s w n s m s u t r s = qpu t r s 
where the product of the space density of sites n s 
and the local site accretion rate m s is set equal 
to qp/tdo and we assume m s /m s = td . Com- 
bining all these results, the effective viscosity for 
the ensemble of sinking cooling sites is roughly 
p s w 1.2 r^® c 5 gm cm -1 s _1 . This viscosity is 
small within most of the bright X-ray core of the 
galaxy. Although p s exceeds p sn for r > 50 kpc, 
there is little or no optical line emission from these 
distant regions. 

Obviously, each of these estimates for p*, p sn 
and p s is highly uncertain, perhaps good to an 
order of magnitude. In addition to these possible 
sources of turbulent viscosity, the interstellar gas 
could become turbulent due to continuous merg- 
ers of gas-rich smaller galaxies having turbulent 
wakes, but there is little optical evidence for suffi- 
ciently numerous merging galaxies. The torquing 
action of magnetic stresses can also propagate an- 
gular momentum outward, masquerading as vis- 
cosity. Finally, we note the considerable energy 
available in differential shear that could be con- 
verted to turbulence if the rotating flow is globally 
unstable. If the characteristic change of azimuthal 
velocity over length is Aw^ « Idv^/dR^^, 
then the effective viscosity is p^ ~ p\dv^j dR^^. 
From Figure 2a, \dv^jdR\ s=s 40 km s _1 kpc -1 
so p<j, ~ lOOOr^ 8 ^ ^, which could be large de- 
pending on the unknown scale length 1$. 



3.2. Viscous Terms 

In rotating cooling flows the flow velocity in 
the meridional plane is very subsonic, typically 
more than an order of magnitude smaller than 
the azimuthal velocity (Paper I). Therefore, only 
the azimuthal velocity component is relevant to 
our investigation of the importance of the viscous 
transport of angular momentum in circularizing 
the X-ray images of elliptical galaxies. In our nu- 
merical solutions, we implement the contribution 
of viscosity by solving the conservation equation 
for angular momentum density ph = pRv^. The 
rate of change of angular momentum density due 
to viscous effects alone is given by 

f d{ph) \ (Ph_ d 2 h p dh 

2h dp dh dp dh dp 
~ ~R ~dR + 'M'dR + ~dzdz 
where p = p(R, z) is a spatially variable viscosity. 
We have incorporated these viscous terms into the 
Eulerian ZEUS2D code (Stone & Norman 1992) 
as a separate computational process. In addition 
to the usual consideration of "source" and "trans- 
port" steps in the Eulerian numerical procedure, 
we modify for viscous transport by solving the 
equation above using a standard operator-splitting 
procedure. The viscosity terms in the angular mo- 
mentum equation above correspond to a modified 
two-dimensional diffusion. We solve these terms 
at each timestep using a second order alternating- 
direction implicit (ADI) difference scheme (Press 
et al. 1992) on the 2D logarithmic grid. Angu- 
lar momentum diffuses spatially because of vis- 
cous effects. When there is no cooling dropout, 
our numerical procedure conserves total angular 
momentum but also allows it to flow out of the 
computational grid at the boundaries. 

The equation for thermal energy conservation 
must also be corrected for the contribution of vis- 
cous dissipation. An additional numerical step is 
applied to the thermal energy equation to solve for 
the following dissipative heating term: 

(^)„„^ fl2 {[4(^)] 2 + [l(^)] 2 } 

where e = 3kT/2p m m p is the specific thermal en- 
ergy. The remaining (non-viscous) terms in the 
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thermal energy equation, and the complete equa- 
tion of motion, are discussed in Brighenti & Math- 
ews (1999a). 

3.3. Effect of Viscosity on X-ray Images 

In general we find that the circularization of X- 
ray images of elliptical galaxies is easier to achieve 
with cooling dropout than with viscous redistri- 
bution of angular momentum. For this reason 
we begin our discussion of viscous rotating cool- 
ing flows with the supernova viscosity, the largest 
viscosity expected from the estimates above. In 
this case \i = H sn {r) where r = [R 2 + z 2 ] 1 / 2 is 
the radial distance from the galactic center. Fig- 
ure 4a shows the X-ray image of our rotating 
galaxy viewed along the equatorial plane at time 
t n = 13 Gyrs but now including viscous terms 
based on /i s „ = 16.2 r^ pc gm cm -1 s _1 . The 
solution shown in Figure 4a also includes a mass 
dropout coefficient q = 1 and the additional sur- 
face brightness T, X (R, z) contribution from the col- 
lective emission from cooling sites. The contours 
in Figure 4a are seen to be almost identical to 
those shown in Figure lc, the same calculation 
without viscosity. The velocity field (not shown) 
is also nearly the same as that in Figure 2b and 
the spatial distribution of dropout gas in this so- 
lution pdo{R,z) is almost identical to the same 
non- viscous q = 1 solution shown in Figure 3. It 
is clear that the additional small viscosity due to 
stellar mass loss and moving cooling sites, /z* and 
/j, s , will have almost no effect on the solutions. We 
conclude that the expected level of turbulent vis- 
cosity is insufficient to circularize the X-ray images 
of cooling flows in elliptical galaxies. 

However, our estimates of the viscosity are ex- 
tremely crude, and it is possible that the real value 
is considerably larger. To explore this idea, we 
show in Figure 4b surface brightness T, X (R, z) con- 
tours at t n = 13 Gyrs computed with a greatly en- 
hanced viscosity, fi — lOOr^' 4 , which still has the 
decreasing radial dependence expected if the vis- 
cosity has a stellar origin. The surface brightness 
contours in Figure 4b, also computed with q = 1 
and dropout emission, are noticeably rounder in 
the bright central region. In addition, a projection 
of the X-ray emission from this solution viewed at 
inclination i = 60° (Figure 4c) has contours that 
are considerably more circular and compares fa- 
vorably with the observations of Hanlan & Breg- 



man (1999). 

Finally, we consider solutions with very large 
and uniform turbulent viscosity. Our estimates 
of the sources of viscosity - fj, sn , /z*, and fi s - 
are independent of galactic shear. However, as 
we speculated earlier, it is possible that these (or 
other) non-linear perturbations could unleash the 
considerable energy density in the shear flow to 
maintain a high level of turbulence throughout the 
flow. Strong, approximately uniform distributed 
turbulence might result from continuing mergers 
of small gas-rich galaxies penetrating through the 
cooling flow. But this hypothesis seems doubtful 
since there is little or no evidence for this type of 
continuous galactic bombardment in either X-ray 
or optical images. If many on-going small mergers 
were required to establish global interstellar tur- 
bulence, numerous small, gas-rich galaxies should 
be apparent in optical images of large elliptical 
galaxies within the optical image and beyond. 

Another obvious motivation for exploring larger 
turbulent viscosities is the persistence of X-ray 
flattening in the 10 - 15 kpc region apparent in 
Figures lb, Id, 4b and 4c, which may still exceed 
X-ray ellipticities in Einstein HRI and ROSAT 
HRI images. Isophotal flattening (and twists) 
have been reported by Boute & Canizares (1998) 
in some large ellipticals, but this effect occurs at 
much larger radii and is probably related to mis- 
alignments of the dark halo and stellar potentials 
which we do not consider here. 

In Figure 5 we show equator-on views of the 
X-ray surface brightness at 13 Gyrs for solutions 
computed with large and spatially uniform viscosi- 
ties, fj, = 5 and 20 gm cm -1 s" 1 respectively. Both 
solutions include dropout with q = 1 and emission 
produced by cooling sites. Clearly, the largest tur- 
bulent viscosity (/z = 20) can circularize images 
both in the 10-15 kpc region and near the bright 
center. Turbulent viscosity has very little influ- 
ence on the current ROSAT luminosity within the 
optical image of the E2 galaxy: L x = 9.1, 11.0, 
10.0, and 11.5 x 10 40 ergs s _1 for q = 1 models with 
the four viscosities we consider, fi = I6.2r^ c 4 , 
100r^" pc , 5, and 10 respectively. 

As we mentioned earlier, the fraction of inter- 
stellar gas that derives from stellar mass loss de- 
creases with galactic radius. At r > 10 kpc a 
larger fraction of the gas has flowed in from the 
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halo (Brighcnti & Mathews 1999a) and is not cre- 
ated solely by stellar mass loss as we have as- 
sumed. The angular momentum density in gas 
flowing in from the outer halo may be different 
(and possibly larger) than that of the stars. How- 
ever, if most of this extended gas was expelled 
from the large elliptical at early times, driven by 
a Type II supernova wind, its angular momentum 
may differ little from that of the galactic stars. 
But, even in this scenario, if this extended gas 
has been tidally influenced by other nearby galax- 
ies, as may be common, when it flows in from the 
outer halo its specific angular momentum could 
differ from the galactic stars both in magnitude 
and direction. 

To explore the possible influence of inflowing 
circumgalactic gas on the apparent X-ray shape of 
galactic cooling flows, we performed several ideal- 
ized calculations in which additional, non-rotating 
circumgalactic gas flows in from the galactic halo. 
In these models (not illustrated here) we used the 
same parameters for the circumgalactic gas as de- 
scribed in Brighenti & Mathews (1998). As ex- 
pected, the ellipticity of the X-ray isophotes at 
10 - 15 kpc is reduced by the inflowing gas, but 
the effect is not very large. Evidently the small 
but significant global spinup resulting from inter- 
action with gas expelled from the rotating stellar 
system is important even at these radii. 

4. VISCOUS DIFFUSION OF METAL- 
LICITY 

If viscosity is invoked to circularize the X-ray 
images of elliptical galaxies or to amplify magnetic 
fields in the interstellar gas (Mathews & Brighenti 
1997), then the level of turbulence required must 
be consistent with observed metallicity gradients 
in the interstellar gas. When X-ray data reduc- 
tion procedures are applied to spatially resolved 
elliptical galaxies, the iron abundance in the in- 
terstellar gas has been found to decrease with in- 
creasing radius approximately as z Fe oc j- -0,5 (e.g., 
Matsushita 1997). We restrict our discussion here 
to iron for which current data is most plentiful. 
In the presence of spatial diffusion the continuity 
equation for iron mass is 



^ + V.(p C u) 



V-(pDVc) 



(1) 



V 1.4 J 



VFeJa 



P 

q^—c, 

tdo 



where c = c Fe = PFe/p is the concentration of iron 
relative to the total mass density. The abundance 
of iron in the hot gas relative to hydrogen mass 
can be found from the concentration by correcting 
for the mass of helium, z = 1.4c. The diffusion 
coefficient depends on a characteristic velocity and 
mixing length, D w vjtd and can be expressed 
directly in terms of the viscosity, D w p/p. Iron 
enters the interstellar gas both from normal stellar 
mass loss and from Type la supernovae. For the 
stellar iron abundance we assume 



z*(R,z) = z Fei0 [l + {R/R„f + {z/z„f 



-0.3 



with central value z Fe , = 1-77 x 10~ 3 (e.g., Ari- 
moto et al. 1997; Ishimaru & Arimoto 1997). 
Each supernova expels yFe.ia = 0.744 Mq in iron 
and a total mass of m sn = 1.4 M Q . With a stel- 
lar mass to light ratio M^/Lb — 9.14, we find 
a sn = 4.85 x 10- 21 SNu(i n )(Vin) _1 s -1 - 

For most cases of interest viscous spatial diffu- 
sion is insufficient to flatten iron abundance gra- 
dients since the gradient is re-established on short 
timescales by Type la supernovae. The timescale 
for spatial diffusion of iron is td, Fe ~ P r2 1 P ~ 
2 x 10 9 r-i p 6 c years using p w 1.05 x lO" 25 ^ 8 gm 
cm~ 3 and p = p sn = lQ.2r^ c gm cm -1 s -1 . By 
comparison, the timescale required for Type la su- 
pernovae to restore an interstellar iron concentra- 
tion of c = z/lA is t sn , Fe « pcm sn /a sn y Fe jaP*- 
Adopting z = 0.5, z Fe ,o = 1-77 x 10~ 3 , p = 



1.05 x lO- 25 ^ 1 / gm cm- 3 , = 4.35 x 10" 



22, 



kpc 



gm cm , and SNu(t„) = 0.03, we find t sn ^Fe ~ 
6 x 10 7 r^, p 2 years. Therefore td 3> t S n,Fe and su- 
pernovae can maintain interstellar iron abundance 
gradients even in the presence of supernova-driven 
turbulence. 

To test this more quantitatively, we have solved 
the iron concentration equation in several evolving 
cooling flows. In cylindrical geometry the diver- 
gence term in the iron continuity equation (1) is 



This diffusion equation was solved using a second- 
order ADI method similar to that used for the vis- 
cosity terms in the equation of motion. 
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Figure 6 compares iron abundance profiles in 
two cooling flows at t = t n along the R and z axes 
for both the gas and stars. For the flow with q = 1, 
P = Msn = I6.2r^ c 4 gm cm -1 s _1 and D = p/p, 
the diffusion has almost no effect on the iron distri- 
bution in the hot gas, in agreement with our sim- 
ple estimate above. Even when the viscosity is in- 
creased to p = 100 throughout the computational 
grid, also shown in Figure 6, an iron abundance 
gradient is still present, but is beginning to flat- 
ten, particularly at the outer edge of the galactic 
stellar distribution at \ogR ~ logz ~ 2. This in- 
dicates that observed interstellar iron abundance 
gradients can be maintained only if the average 
turbulent viscosity is p < 100 gm cm -1 s _1 . To 
investigate further the importance of the source 
terms on the metallicity gradient, we performed an 
additional exploratory computation (with q = 1 
and p = 100) in which the iron source terms in 
equation (1) were set to zero beginning at time 
t = 10 Gyrs. By time t = 11 Gyrs in this calcu- 
lation the iron gradient in the interstellar gas was 
completely leveled, resulting in uniform ZFe(R, z) 
throughout the computational grid. In this unre- 
alistic test calculation our solution of the concen- 
tration diffusion term above behaved as expected. 

The diffusive mixing of supernova enrichment 
in the interstellar gas is relevant to the possibil- 
ity that locally metal-rich regions may cool by X- 
radiation and dropout faster than the rest of the 
cooling flow. We estimate the metallicity diffusion 
time by equating the diffusion length from each su- 
pernova blast wave bubble to the mean distance 
between bubbles. The number of Type la super- 
nova bubbles in a large elliptical galaxy at any 
time, Af bubb s=s v sn t huhh ~ 1.5 x 10 3 (t 6ttb6 /10 6 yrs), 
is quite small. Here tbubb ~ 10 6 yrs is a typ- 
ical lifetime of a supernova bubble to destruc- 
tion by buoyant mixing (Mathews 1990). How- 
ever, the iron deposited in the cooling flow by 
Type la supernova persists long after the ther- 
mal energy in the hot blast wave bubble has dis- 
sipated. To estimate the efficiency that super- 
nova iron mixes throughout the interstellar gas 
by turbulent diffusion, we imagine a time t = 
at which the supernova explosions begin. Then 
the mixing time t m i X for supernova iron to dif- 
fuse between supernova sites is found by equating 
the decreasing distance between supernova bub- 
bles at t > to the distance that iron has dif- 



fused by turbulent mixing from any particular site; 
tmix must then be compared to the local cooling 
flow time. As Type la supernovae events increase 
with time t, the mean separation between sites is 
hubb ~ {tv S nP*/M itt )~ 1 / 3 ~ 36r kpc tyr /3 kpc, esti- 
mated with p* = 4.35 x 10~ 22 r^ pc gm cm~ 3 . This 
can be compared with the distance Sdi / / that dif- 
fusion mixes metals from an individual supernova 
bubble in the same time t: Sdiff — (pt/p) 1 ^ 2 ~ 
2.3 x ^0~ 5 tyr r kpc kp c where we assume p = p sn = 
I6.2r-^ gm cm" 1 s" 1 and p = 1.05 x lO" 25 ^ 8 
gm cm~ 3 . Diffusive mixing of supernova enrich- 
ment in the interstellar medium occurs in time 
t = Unix found by setting 5 dlf f « S bu bb, which 
gives t mix ~ 3 x 10 7 r° p 9 c 6 yrs. This mixing time 
is comparable to the radial flow time in the cool- 
ing flow, tfiow ~r/«~5x 10 7 r fe p C [u/20(km/s)] _1 
yrs. We conclude that supernova enrichment prod- 
ucts are well mixed throughout the cooling flow for 
the Type la supernova rate we assume here; for 
turbulent viscosities p > p sn , mixing is guaran- 
teed. In a cooling flow the local radiative cooling 
time is comparable to the flow time, so iron pro- 
duced by supernovae disappears from the flow by 
mass dropout in time ~ tfi ow . The complete ab- 
sence of diffusive mixing, as envisioned by Fujita, 
Fukumoto & Okoshi (1996; 1997), in which un- 
mixed, metal-rich supernova bubbles and stellar 
ejecta cool rapidly by enhanced X-ray line emis- 
sion and differentially drop out from the flow may 
not be realistic, but we have not explicitly ruled 
out this possibility 

5. SUMMARY AND CONCLUSIONS 

In this paper we have addressed the curious ab- 
sence of X-ray disks in rotating, luminous ellipti- 
cal galaxies. In our previous hydrodynamic stud- 
ies of rotating cooling flows (Paper I; Brighenti 
& Mathews 1997) both angular momentum and 
mass were conserved and the resulting X-ray im- 
ages were dramatically flattened out to 1-2 effec- 
tive radii when viewed perpendicular to the axis 
of rotation. These results appear to be inconsis- 
tent with the more circular X-ray images found by 
observations with Einstein HRI (Fabbiano, Kim 
& Trinchieri 1992) and ROSAT HRI (Hanlan & 
Bregman 1999). Nevertheless, the X-ray shapes of 
elliptical galaxies are difficult to determine from 
the observations and we may need to wait for the 
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next generation of X-ray telescopes with higher 
spatial resolution for definitive determinations of 
X-ray ellipticities. It is likely that interstellar ro- 
tational flattening will eventually be observed. 

In the preceding discussion we have explored 
the circularizing influence of mass dropout and 
turbulent transfer of angular momentum on the 
X-ray images of rotating galactic cooling flows. 
Strong theoretical and observational arguments 
support both of these possibilities. Localized ra- 
diative cooling leading to mass dropout through- 
out the cooling flow is indicated by the relatively 
low masses of central black holes and adjacent nu- 
clear regions when compared to the mass of in- 
terstellar gas that is expected to cool over cos- 
mic time. Further evidence for cooling dropout is 
provided by optical line emission within R e and 
the evidence in this same region for a young stel- 
lar population having masses that extend up to 
1-2 M Q but not beyond (Ferland, Fabian, & 
Johnstone 1994; Mathews & Brighenti 1999a). In- 
terstellar turbulence is also very likely in ellipti- 
cal galaxies, generated by mass transport associ- 
ated with stellar mass loss, Type la supernovae 
and mass dropout. Indirect evidence for cooling 
flow turbulence is provided by interstellar mag- 
netic fields of several fiG at <~ 10i? e , which are 
typically observed in elliptical galaxies having dou- 
ble radio sources. Fields of this magnitude can be 
understood as originating from small seed fields 
ejected from mass-losing stars, followed by subse- 
quent field amplification by an interstellar turbu- 
lent dynamo mechanism (Moss & Shukurov 1996; 
Mathews & Brighenti 1997). Due to the field con- 
centration associated with the inward, converging 
motion of the cooling flow, interstellar magnetic 
fields are further magnified and are expected to be 
particularly strong, perhaps exceeding equiparti- 
tion values, in the central regions of cooling flows. 
Owen and Eilek (1998) find fields of 10 - 100 ^G 
within r = 50 kpc in the bright cD galaxy NGC 
6166. 

Elliptical galaxies of lower optical luminosity 
Lb < 3 x 10 10 Lb.q differ in many qualitative 
ways from the more massive elliptical galaxies sim- 
ilar to the one we have studied here (Faber et 
al. 1997). Low luminosity elliptical galaxies ro- 
tate faster and are rotationally flattened. We have 
shown (Brighenti & Mathews 1997) that rotating, 
non-dropout cooling flows in these galaxies also 



form massive, extended disks of cooled gas. How- 
ever, when q = 1 mass dropout is included in these 
calculations, disks of cold HI (or H 2 ) do not form 
from the rotating cooling flow. There is observa- 
tional evidence for extended disks of cold HI gas in 
some low luminosity elliptical galaxies, but these 
HI disks often extend far beyond the optical im- 
ages of the galaxies, indicating that they were cre- 
ated by a (possibly very old) merging event rather 
than by cooling flow dropout (Oosterloo, Morganti 
& Sadler 1999a; 1999b). In some cases the cold 
gaseous disks have central holes in neutral hydro- 
gen replaced with HII emission (Oosterloo, Mor- 
ganti & Sadler 1999), suggesting a more complete 
conversion into stars in this part of the disk where 
gas pressures and densities are highest and star 
formation (and ionization by PAGB stellar UV) 
should be most efficient. In general, systematically 
larger H/3 features are observed in the stellar spec- 
tra of elliptical galaxies with Lb <3x 10 10 Lb,©, 
indicating considerable recent star formation (de 
Jong & Davies 1997). Unfortunately, thermal X- 
ray emission from the hot interstellar gas in these 
low-luminosity, low-E^ ellipticals is masked by the 
collective emission from low mass X-ray binaries. 
A detection of X-ray rotational flattening in their 
cooling flows will require spectral separation of the 
hot gas and stellar X-rays with high spatial resolu- 
tion, clearly a job for the next generation of X-ray 
telescopes. 

When mass dropout and turbulence arc con- 
sidered in rotating cooling flows inside luminous 
elliptical galaxies, we find that: 

(1) Mass dropout alone strongly circularizes the 
X-ray images in rotating cooling flows since 
gas is removed from the flow before it has 
moved very far from its point of origin (stel- 
lar mass loss) toward the axis of rotation. 

(2) Conversely, the absence of strong rotational 
flattening in X-ray images of elliptical galax- 
ies is persuasive evidence for distributed 
mass dropout. 

(3) In the presence of mass dropout by radiative 
cooling with q = 1, no cold gaseous disks 
form on the equatorial plane having radii 
larger than our innermost grid size, ~ 150 
pc. This is due to the sensitivity of radiative 
cooling (and associated mass dropout) to the 
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local density, {dp/dt)d oc qp/tdo oc p 2 . 



(4) The spatial distribution of cooled gas mass is 
markedly natter than that of the old stellar 
population, especially near the galactic cen- 
ter. Any optical signature of a younger stel- 
lar population formed from the cooled gas, 
such as the H/3 index, should also exhibit 
a significantly higher ellipticity than that of 
the older, background stars. 

(5) The estimated viscosity from known sources 
of interstellar turbulence is dominated by 
motions induced by Type la supernovae; 
stellar mass loss and mass dropout are 
smaller sources of interstellar turbulent vis- 
cosity. 

(6) Supernova-induced turbulent viscosity is in- 
sufficient to circularize the X-ray appearance 
of rotating elliptical galaxies; the viscous ef- 
fects of stellar mass loss and cooling dropout 
on the X-ray images are even smaller. 

(7) Much larger, spatially uniform turbulent vis- 
cosities can circularize the X-ray isophotes 
throughout the galaxy; such turbulent vis- 
cosities could result from rotational shear 
instabilities, but we have not demonstrated 
this instability here. 

(8) The spatial diffusion of interstellar iron due 
to turbulent viscosity docs not appreciably 
reduce the interstellar iron abundance gra- 
dients observed in bright elliptical galax- 
ies, provided the viscosity does not exceed 
~ 100 gm cm -1 s _1 . This follows since the 
observed interstellar iron gradients are re- 
established on short time scales by Type la 
supernovae and, to a lesser extent, by stellar 
mass loss. 
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Fig. 1. — X-ray surface brightness distribu- 
tions T, x (R,z) (erg cm~ 2 s _1 ) for emission in 
the ROSAT band (0.1 - 2.4 keV) shown in one 
quadrant of the rotating E2 galaxy at t = 13 
Gyrs. Several contours are labeled with values 
of log£ x (ii, z). (a): Rotating cooling flow with 
no mass dropout (q = 0) viewed perpendicular to 
the axis of rotation; contours are evenly spaced by 
A(log£ x ) = 0.15. (b): Same cooling flow shown 
in (a) but viewed at an inclination of 60°. (c): 
Rotating cooling flow with mass dropout (q = 1) 
viewed perpendicular to the axis of rotation; con- 
tours are evenly spaced by A(log£ x ) = 0.3. (d): 
Same cooling flow shown in (c) but viewed at an 
inclination of 60°. 

Fig. 2. — Azimuthal velocity distribution v ( j ) (R, z) 
(km s _1 ) in the cooling flow at t = 13 Gyrs when 
viewed perpendicular to the axis of rotation, (a): 
v^(R, z) for a cooling flow with no mass dropout 
(q = 0); the spacing between adjacent contours is 
29.4 km s _1 , which is also the velocity on the con- 
tour closest to the z-axis. (b): v ( f,(R, z) for a cool- 
ing flow with mass dropout (q = 1); the spacing 
between adjacent contours is 23.5 km s _1 , which 
is also the velocity on the contour closest to the z- 
axis. (c): the emission- weighted line of sight veloc- 
ity vi os (R 7 z) for a cooling flow with mass dropout 
(q = 1); the spacing between adjacent contours 
is 14.7 km s _1 , which is also the velocity on the 
contour closest to the z-axis. 



T, x (R, z) computed with viscosity fi sn = 100r fc ' 
gm cm -1 s _1 and viewed perpendicular to the axis 
of rotation, (c): T, X (R, z) computed with viscos- 
ity H sn = lOOr^' 4 gm cm -1 s _1 but viewed at 
inclination i = 60°. 

Fig. 5. — ROSAT band X-ray surface brightness 
distributions T, X (R, z) (erg cm -2 s _1 ) in the rotat- 
ing E2 galaxy viewed perpendicular to the axis of 
rotation at t — 13 Gyrs. Several contours are la- 
beled with values of log T, x (R, z) and the contours 
are evenly spaced by A(logS a: ) = 0.3. (a): Com- 
puted with uniform viscosity \l — 5 gm cm -1 s _1 . 
(b): Computed with uniform viscosity \i = 20 gm 
cm -1 s _1 . 

Fig. 6. — Interstellar iron abundance profiles along 
the R axis (solid lines) and the z axis (dashed 
lines) at 13 Gyrs for the rotating E2 galaxy with 
q = 1 and two viscosities: /U = /i s „ = 16.2r^ p 1 c 4 
(heavy lines) and uniform fi = 100 (light lines). 
The dotted lines show the stellar iron abundance 
along the R axis (heavy dotted line) and the z axis 
(light dotted line). 



Fig. 3. — Comparison of the spatial distribution 
of cooled gas (and young star) density pdo(R,z) 
(solid contours) with that of the old stellar pop- 
ulation p*(R,z) (dashed contours) for the q = 
1 solution at 13 Gyrs. Contours are labeled 
with \ogpdo(R, z) or logp*(i?, z) with units gm 
cm~ 3 . The spacing between adjacent contours is 
A(logp) = 0.35. 

Fig. 4. — X-ray surface brightness images T, X (R, z) 
of the E2 galaxy at t n = 13 Gyrs in the ROSAT 
band including emission from mass dropout with 
q = 1. Several isophotal contours are labeled with 
log Yi x in cgs units and contours are evenly spaced 
in A(logSa;) = 0.3. (a): T, X (R, z) computed 
with viscosity fi sn — 16.2r^ c 4 gm cm -1 s _1 and 
viewed perpendicular to the axis of rotation, (b): 
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Table 1 

GALACTIC PARAMETERS FOR 
EO GALAXY 



Parameter 


Value 




311.59 pc 




5.088 kpc 


Rt 


113.1 kpc 


ft = Rch/Rc* 


19.953 




1.438 x KT 20 gm cm" 3 


7 = Pho/p*o 


3.7954 x KT 4 


M H 


4.52 x 10 11 M 


M ht 


9 


L B 


4.95 x 1O 1O L B0 


M* t /L B 


9.14 


cr» b 


351 km s" 1 



a Effective radius. 

b Characteristic velocity dispersion in 
stellar core, cr* = (47rG/?* i? 2 < ,/9) 1 ' /2 . 
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